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Abstract 

Following the nonrelativistic QCD factorization scheme, by taking latest available measure¬ 
ment on Xbi^P) into consideration, we present an updated study on the yield and polarization 
of T(1S', 25*, SS") hadroproduction, and the fractions of Xbi'bnP) feed-down in T{nS) production 
at QCD next-to-leading order. In the htting, three schemes are applied with different choice of 
XbiiTiP) feed-down ratios and NRQCD factorization scale. The results can explain the measure¬ 
ments on yield very well as in our previous work. The polarization puzzle to T (35) is now solved by 
considering the Xb{3P) feed-down contributions. The ratio of a\xb 2 {lP)\/cr[xbi{lP)] measurements 
from CMS can also be reproduced in our prediction. Among the different schemes, the results show 
little difference, but there are sizeable difference for the fitted long-distance color-octet matrix el¬ 
ements. It may bring large uncertainty when the values are applied in other experiments such as 
in ee, ep colliders. 

PACS numbers: 12.38.Bx, 13.60.Le, 13.88.-he, 14.40.Pq 
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I. INTRODUCTION 


Quantum chromodynamics(QCD) successfully describes the strong interaction at parton 
level due to its asymptotic freedom property. But it fails to calculate observations with 
detected hadrons directly since the hadronization from quark is nonperturbative. There¬ 
fore, factorization scheme to bridge the perturbative calculable part and nonperturbative 
hadronization part is crucial. For heavy quarkonium production and decay, non-relativistic 
QCD(NRQCD) factorization schemejl, 2], which was proposed to explain the huge discrep¬ 
ancy between the theoretical prediction and experimental measurement on the transverse 
momentum distribution of J/'ijj production at the Tevatron, has been a very successful 
scheme in many applications. However, it encounter challenges on the transverse momen¬ 
tum distribution of polarization for J /ip and T hadropro duct ion where the theoretical pre¬ 
dictions can not describe the experimental measurements at QCD leading order (LO), or 
next-to-leading order (NLO). 

In the last few years, signihcant progresses have been made in NLO QCD calculat 
The NLO corrections to color-singlet J/ip hadropro duct ion was investigated in Refs. 


ions. 




where the pt distribution is found to be enhanced by 2-3 orders of magnitude at the high pt re¬ 
gion and the J/tp polarization changes from transverse into longitudinal at NLO jd]. And the 
results are reproduced at LO in a new factorization scheme for large pt quarkonium produc¬ 


tion 


j^. In Ref. jo], NLO corrections to J/xp via S-wave color-octet (CO) states(^S'Q^^, ^S'f^)are 


studied and the pt distributions of both J/%p yield and polarization changec 
with LO. The NLO corrections for Xcj hadropro duct ion are studied in Ref. 


little compared 
7|. Calculations 


and hts for bot 
three groups 


i_j^ield and polarization oi J/tp production at NLO QCD are presented by 
Where the complete prompt J/pj hadroproduction study, which is cor¬ 
responding to the present available experimental measurements, is performed in Ref. for 


the hrst time. Th e J/ip polarization puzzle is still not fully understood. Recently, the LHCb 

-141 came 


Collaboration [u| published their measurement of rjc production. Three letters 
out successively, investigating the data from different points of views. Ref. [l2| considered 
the Pc experiment as a challenge of NRQCD, while Refs. [l3, 14| emphasized its indications 
on the J/pj productions and polarizations. The complicated situation suggests that, further 
study and phenomenological test of NRQCD is still an urgent task. 

T production and polarization is an alternative one of the best laboratories for under- 
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standing the physics in the hadronization of heavy qnark pairs. Dne to its heavier mass 
and smaller v (where v is the velocity of the heavy qnark in the meson rest frame), one can 
expect better convergence in the QCD and NRQCD expansions, and conseqnently better 
description of the experiment by QCD NLO predictions. For T hadroprodnction, similar 
progresses are also achieved on the pt distribntion of y ield and polarization for the CS chan¬ 


nel at QCD NLO 3,_4|, and for S-wave CO states 15|. The NLO correction via all CO state 


is stndied in Ref. [1^. The hrst complete NLO QCD corrections on yield and polarization 
of T(1S', 2S', SR) are presented in onr former work 17|, where the resnlts can explain the 


experiments on yield very well, so as the polarization of T(1S', 2S') at CMS measnrements. 
However, withont considering the Xbi^P) feed-down contribntion, the polarization of T(3S') 
is inconsistent with data. Thereafter, the mass of xti^P) is measnred 18j at the LHC, and 
the fractions for T(3S') prodnction from Xb(3P) radiative decay is hrst measnred [l9| by the 
LHCb Collaboration. The large measnred valne of this fraction indicates that the reconsid¬ 
eration for T(3S'), as well as T(1S', 2S'), is needed. In this work, we take into consideration 
of Xfe(3P) feed-down contribntions carefnlly to npdate the yield and polarization analysis on 
T(nS') hadroprodnction at QCD NLO correction, and also predict the ratio of differential 
cross sections of Xb 2 i^P) to that of Xbii^P) on the LHCb 18| and CMS data 20|. 

This paper is organized as follows. In Sec. im we give a brief description of the theory in 
onr work. Nnmerical resnlts are presented in Sec. IHII In Sec. IIVI we introdnce more htting 
schemes and discnss the difference between them. Finally, the snmmary and discnssions are 
given in the last Section. 


II. THEORY DESCRIPTION 

Following the NRQCD factorization formalism j^, the cross section for qnarkoninm 
hadroprodnction H can be expressed as 


da\pp H + X] = / dxidx 2 GpG^p x dd[ij {bb)nX]{0^) 


( 1 ) 




where p is either a proton or an antiproton, is the parton distribntion fnnction(PDF) 
of p, the indices i,j rnns over all possible partonic spices, and n denotes the bb intermediate 
states (^S'fVS'o®',^S'f',^Pj®') for T and (^Pj^',^S'f') for XbJ- The short-distance contribntion 
dd can be calcnlated pertnrbatively, while the long-distance matrix elements(LDMEs) (O^^) 
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are governed by nonperturbative QCD effects. 


The polarizations of T are defined as 


2l| 


dan — dcToo \/2Redaio y/2Redai _i 

A = ^^^—, u = 


( 2 ) 


dan T da^Q dan T da^Q dan T dcoo 
where das^si is the spin density matrix of T hadropro duct ion. Only the parameter A in 
helicity frame is considered in our work. 

The fractions of T{nS) originating from Xb(mP) decays is dehned as 


n 


Xb{mp) _ 
T(n5) — 




a{pp Xbj{mP)X) 
a{pp —)■ T(nS')X) 


X B[xb{mP) T(n5)] 


( 3 ) 


J=0,l,2 

where n and m are radial quantum numbers of the bound states and B denotes the branching 
ratios for the decay Xbj{^P) —t T(nS') 7 . 

To obtain das^s'^i similar treatment as in Ref. [w, 17| is taken for both direct and feed- 
down contributions. For various feed-down contributions in T production, we treat them in 
different ways to get comparable results. The details are given in the next section. 

The newly updated Feynman Diagram Calculation package 2^, [2^ and the are used in 
our calculation. Compared with our former work, the newly added calculations are about 
the productions on Xb{^P) in different experimental conditions. 


III. NUMERICAL RESULTS 


In the numerical calculation, the CTEQ6M parton distribution functions 2^ and cor¬ 
responding two-loop QCD coupling constants as are used. We adopt an approximation 
rub = Mh/2 for the quark mass, where Mh is the mass of bottomonium H. All the masses 


are taken from PDG 


25[ |. except for Xbj{^P), which are chosen as =10.511 GeV 


for J=0,l,2 19|. Therefore, the mass of bottom quark in our calculation is different values 
for T(nS') and XbipP)- 

The color-singlet LDMEs are estimated from wave functions at the origin 




( 4 ) 


while the wave functions and their derives via potential model calculation 2^. We listed 
the results in Table HI 
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TABLE I: Radial wave functions at the origin 


T{nS) 

PT(nS)(0)P 

XbimP) 

RA&(mn)(0)P 

IS 

6.477GeF3 

IP 

lAllGeV^ 

2S 

3.234GeF3 

2P 

1.653GeF^ 

3S 

2.474GeF3 

3P 

L794GeF'^ 


Branching ratios involving bottomonia can be found in Ref. [17[ Table I, in which however, 
that for Xbj(3P) are not included. Since no experimental data for branching ratios of Xfe(3P) 
feed-down to T(?7,S') is available right now, we take B[xbj{^P) T(3S')] ~ B[xbj{‘^P) —t 
T(2S')] as an approximation and ignored the contributions from Xbj{^P) for T(2S') and 
T(IS') due to the little fractions. 

The factorization, renormalization and NRQCD scales are chosen as pj = = 

and /ia = fn^v ~ 1.5 GeV, respectively. The center-of-mass energy is 1.8 TeV 
and 1.96 TeV for Tevatron Run I and Run II and, 7 TeV and 8 TeV for the LHC, respec¬ 
tively. A shift ~ p^' X {Mh/Mh>) is used while considering the kinematics effect in the 
feed-down from higher excited states. 


In the fit, we have used three kinds o. 
the differential cross section from CDF 


data from experimental measurements, namely. 


P, LHCb |28|, CMS |29| and ATLAS |30|; the 


polarization from CDF 3]| and CMS |32|; the fractions of T(nS') production originating 
from radiative decays of x;,(mP) meson (Px*(n 5 ^^) Ql- included the data of the fraction 
from Xb(3P) feed-down at both ^/s=7 TeV and \/i=8 TeV in the fit of T(3S'), while for 
T(2S') and T(IS'), only the points at \/i=7 TeV for fractions are included. The linear 
interpolation method has been taken when dealing with the fraction data since it behaves 
smoothly and the pt points nearly our theoretical calculations are chosen. Only the data in 
the region pt > 8 GeV are used in our fit as we know that the double expansion in as and 

is not good in the small pt regions. 

We perform the fits for T(3S'),T(2S'),T(1S') hadroproduction step by step, and the cor¬ 
responding x^/d-o.f are 97172^ 114/47, 73/44. All the fitted CO LDMEs are presented in 
Table m A covariant-matrix method [1^ is performed for the plots in order to express the 
uncertainty from the CO LDMEs properly. But we only rotate the direct three LDMEs, 
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TABLE II: The obtained CO LDMEs for bottomonia production(in units of 10“^ GeV^). 
Xbi^P) feed-down contributions are only considered for T(3S') with branching ratios B[x6j(3P) —)• 
T(3S')] ~ B[xbj{‘^P) —^ T(25)] as an approximation. 


state 

(C>T("S)(l5'[8])) 

(C>^(”^)(3S[®1)) 


state 

(^QXbimP) ^3 

Tils') 

13.6 ± 2.43 

0.61 ± 0.24 

-0.93 ± 0.5 

XbilP) 

0.94 ± 0.06 

T(2S) 

0.62 ± 1.98 

2.22 ± 0.24 

-0.13 ± 0.43 

Xbi2P) 

1.09 ± 0.14 

T(3S) 

1.45 ± 1.16 

1.32 ± 0.20 

-0.27 ± 0.25 

Xbi^P) 

0.69 ± 0.14 


namely, while the last one Sf)) is 

fixed. The descriptions for yield, polarization and fractions are given in the following sub¬ 
sections. 

A. Yield and polarization 

The results for differential cross section of T hadroproduction are shown in Fig. [1], while 
those for polarization are shown in Fig. |2l The uncertainty bands in the hgures are from the 
error of CO LDMEs. The experimental data are stamped with different energy y/s, rapidity 
cut and corresponding collaboration in all pictures. 

We can see from the Fig. [1] that the results on the yield of T(nS') hadroproduction £t all 
the experimental measurements very well in a wide pt range and the uncertainties are very 
small. In contrast with our previous £t results jl^ without Xb{^P) feed-down contributions, 
where the polarization of T(3S') shows a weird steep behaviour and can not explain the 
experiment measurement. Our updated hts for T(3S') has changed a lot with the contribu¬ 
tion of Xfe(3P) feed-down for the polarization shown in Fig. |2l and it hts the experimental 
measurements very well. For T(2S'), the production is dominated by channel, which 
lead to a slightly transverse behaviour(A ~ 0.6) at high pt region. Nevertheless, the results 
for polarization can explain the CMS data well, but the distance from the CDF experimen¬ 
tal data is still there. This is also the case for T(IS') polarization, although it is almost 
unpolarized in all pt region (0 < A < 0.4). 
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FIG. 1: Differential cross section for T hadroproduction at the Tevatron and LHC. From left to 
right: T(15), T{2S) and T{3S). Rows from top to bottom correspond to different experimental 



B. Fractions 

In Fig. [31 we give the results on the fraction of Xhi'f^P) feed-down to T(nS'). Only 
four fractions ^ ^ presented since in our approximation the 

Xbi^P) feed-down contributions are ignored for T(1S', 25') production. For , we give 

the results at ■s/s=7 TeV (black solid line) and s/s =8 TeV (red dotted line), but they are 

almost overlapping. For the other three cases, the results at Ps =8 TeV are not included 

□ 

since in the experimental measurements they are almost the same with the ^/s =7 TeV [19[ 
and it is believed that theoretical predictions for the fractions is almost unchanged for the 
centre-of-mass energy ^/s =7,8 TeV. Our results for show different behaviours for 

each n = 1,2,3, they are increasing, flat {pt >10 GeV), slightly decreasing respectively as 
Pi increases. And all of them £t the experimental data very well. The remaining fraction 
shows a increasing behaviour but the theoretical result overshoot the data about 
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FIG. 2: Polarization parameter A for T hadroproduction at the Tevatron and LHC. Prom left to 
right: T(1S), T(25) and T(3S'). Rows from top to bottom correspond to different experimental 
conditions of CDF run II, CMS(|?/| < 0.6), and CMS(0.6< \y\ <1.2). The experimental data are 


taken from Refs. 
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FIG. 3: Fractions functions oi pj . From left to right: ™ 

and ^xT 5 )^’^TTS)^ ™ second row. The experimental data are collected from Ref. {l^ . 
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C. Ratio of cr(xfe2(l-P))/o-(X6i(l-P)) 


Recently, new measurements on the ration of cross sections of a[xb 2 {^P)] to (T[xbi(lP)] 


are reported by LHCb 


181 and CMS collaborations 201. 



FIG. 4: The ratio of the cross sections of Xb 2 {^P) to Xbii^P) production , as a function of pt- 
The yellow band is the results for CMS experimental condition while the green band is for LHCb 
experimental condition. The data are taken from Ref. 


18| 




As a way to check our htting results for the LDMEs, we present the prediction of the 
cross section ratio cr[xb 2 (lP)]/o'[xfci(lP)] in Fig. IH where the green band is the results for 
LHCb experimental condition while the yellow band is for CMS condition. As the picture 
shows that the CMS experimental data can be explained by the theory, while the LHCb data 
is underestimated. One may notice that here is an approximation for the CMS theoretical 
prediction that the rapidity range is \y\ < 1.2, while for the experimental data in Ref. 
the rapidity range is \y\ < 1.5. Although the prediction condition is not exactly the same 
with data, we think the approximation is fairly enough. 


IV. MORE FITTING RESULTS 


While we are preparing this work, a similar work was preprinted by H.Han et al. 33|, 


where they estimated the branching ratios of Xbi^P) feed-down to T(1S', 2S, 35) and htted 
the productions and polarizations. In our £t in the former section, we ignored the Xbi^P) 
feed-down contributions for T(15, 25) since it is really small. To give a more rigorous feed- 
down contributions would be instructive. Meanwhile, the NRQCD factorization scale 
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dependence dose exist at the fixed order results. To explore the uncertainty from these two 
aspects, we refit the production and polarization and obtain two set of new LDMEs for 
Pa = TnfjV = 1.5 GeV and pa = together with the branching ratios B[xbj{3P) —^ ^(^>5')] 


estimated in Ref. 


33| TABLE I, and the numerical results are presented in Table. IIIII and 


Table. HVl respectively. 

For convenience, we define three short statements. 


default scheme (pa = rribV = 1.5 GeV and our naive estimation for branching ratios 
B[xb{^P) ^ T(3^)] = B[xb{2P) ^ T(2^)], B[xb{^P) ^ T(l^,2^)] = 0), 

Han’s scheme (pa = mi, and branching ratios B[xbj{^P) ^(’^>5')] estimated in 


Ref. 


33| TABLE I), 


mix scheme (pa = mhV =1.5 GeV and branching ratios B[xbj{3P) —^ esti¬ 

mated in Ref. [331 TABLE I). 


For the mix scheme in Table. IIIII some changes compared with default scheme in Table. ITT] 
should be mentioned. For T(3S'), the dominate GO channel becomes and the numerical 
value for becomes larger. This is not at all a surprise since we know the 

branching ratios B[xb{^P) T(3S')] used in the default scheme is much larger than it is in 
the mix scheme. This difference is consistent with the change of the proportions for LDMEs 
in two cases. For T(2S') and T(IS'), the LDMEs values have smaller changes. 

In comparison with the mix scheme, for Han’s scheme in Table. IIVI the LDMEs for T(2S') 
changed a lot because the channel becomes negative and the ^Pf^ channel gives positive 
contributions. Besides, the value of becomes larger, which may impact the 

behaviour of the ratio a[xb 2 {)-P)]/cr[xbi{)-P)]- 


TAB 

Ref. 


E III: Same as Table. m except that we choose the branching ratio B[xbj{^P) T(nS')] in 
33|, so-called mix scheme. 


state 


(DT(™S)(35.[8])^ 

(D^("^)(3R|®]))/m2 

state 

(^QXb(mP) (3 

T(IS') 

10.1 ± 2.23 

0.73 ± 0.22 

-0.23 ± 0.50 

Xb{lP) 

0.91 ± 0.06 

T(2S') 

1.19 ± 1.93 

1.88 ± 0.23 

-0.01 ± 0.42 

Xb{2P) 

1.07 ± 0.12 

T(3S') 

-0.15 ± 0.90 

1.53 ± 0.12 

-0.02 ± 0.19 

Xbi^P) 

1.76 ± 0.14 
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TABLE IV: Same as Table. IIIIl except that we choose //a = ^6) the Han’s scheme. 


state 


(0T(nS)(3^[8])^ 

(D^("D(3p[8]))/^2 

state 


T(IS') 

11.6 ± 2.61 

0.47 ± 0.41 

-0.49 ± 0.59 

XbilP) 

1.16 ± 0.07 

T(2S') 

-0.59 ± 2.31 

2.94 ± 0.40 

0.28 ± 0.52 

Xb{2P) 

1.50 ± 0.21 

T{3S) 

-0.18 ± 1.40 

1.52 ± 0.33 

-0.01 ± 0.30 

xR3P) 

1.92 ± 0.34 


The results of the yield and polarization in both Han’s scheme and the mix scheme are 
almost the same as those in default scheme, and the yield results describe the experimental 
data very well and so as the polarization of CMS measurement, but not the CDF measure¬ 
ments. We would not show these pictures here to avoid unnecessary repetition since there 
is little difference compared with results in default scheme shown in Fig. [1] (for yield) and 
Fig. [2] (for polarization). 

In Fig. [5l we presents all the results for the fractions of XbijnP) feed-down to T{nS) 
obtained by using the LDMEs in the mix scheme and Han’s scheme. The green bands are 
our htted results while the yellow bands are obtained by using the LDMEs in Ref. j^. 

We summary the results in two parts. The First part, our results show little difference 
for fractions between Han’s scheme and the mix scheme, except for which in Han’s 

scheme the band is more close to the data. Since the only difference between the two fitting 
scheme is the choose of NRQCD factorization scale /ta, it indicates the scale dependent for /xa 
is small to the hxed order correction. The Second part, we show the difference between our 
results and H.Han’s results in Ref. [33|. In their results, the XbijnP) feed-down contributions 
tend to increase as pt goes higher for all T(?7,S'). In our results, however, only the fractions 
keep increasing behaviour for all the three schemes, and the other fractions 
show a slightly decreasing or smooth behaviour. Meanwhile, with a large uncertainty band, 
H.Han’s results can cover all the experimental data except for Xb(lP) to T(IS'), where their 
results underestimated the data gently. Our results present a somewhat interesting case, 
that the results in our three schemes can give a good explanation to the fractions 
for m = n=l,2,3, while overestimated the data by a factor less than 2 for the other fractions 
that m ^ n. 

We present the ratio of cross section a[xh 2 i)-P)\l(^[Xbii)-P)] for Han’s scheme in Fig. [6l 
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FIG. 5: Fractions functions of The green bands are our predictions that, the upper 

two rows are the results for Table. m mix scheme while the lower two rows are for Table. HVl Han’s 
scheme. The yellow bands are obtained by using the LDMEs in Ref. 3^. From left to right: 

™ third row and ™ second 

R 


"'T(3S) ’ '^T{2S) ’ '^T{2S) tiiiiu luw ain.i 

and forth row. The experimental data are collected from Ref. 


while the results for the mix scheme are not presented to avoid repetition since they are 
indistinguishable with that for default scheme shown in Fig. IHdue to the almost the same 
value of (see Table. UTland Table. [ITT|) . 

It can be seen that the value of (C)^*’T-P)(35'[®]^^ Table. [IV] for Han’s scheme is larger 
than that in Table. HI for default scheme. With this value, the Xb 2 {^P) and x;,i(lP) ratio 
predictions in Fig. |6]for Han’s scheme describe the experimental measurements better than 
that in Fig. 0] for default scheme, and the results cover all the CMS data very well while 
only one point for LHCb measurements is inside the error band. Nevertheless, this difference 
can be explained as the uncertainty of NRQCD factorization scale. In our results for both 
CMS and LHCb, all the centre value of the ratio a[xh 2 {XP)]/cr[Xbi{)-P)] are less than 1, or 
equally, there is a relation cr[xfe 2 (lP)] < <^[Xfei(l-P)]. This is the case for CMS experimental 
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FIG. 6: The same as Fig. HJ except that the LDME in Table. HV] for Han’s scheme 

is used . 


data. But for LHCb, the value at two points of three are bigger than 1. In this point of 
view, the prediction might be in inconsistent with LHCb data. However, we can not make 
the conclusion for the ratio here due to the few LHCb data. 


V. SUMMARY AND CONCLUSION 


y on the yield and polarization of 
171 without Xbi^P) feed-down contri- 


In this work, we presented an updated stuc 
T{1S,2S,3S). In contrast to our previous study 
bution, there are two new points: Xb{^P) feed-down contribution is taken into consideration 
by using available measurement on Xb{3P)] the experimental measurements on the fraction 
of Xb(2P, IP) feed-down is available and applied in our £t. We obtain the CO LDMEs for 
T hadroproduction by htting the experimental data for yield, polarization and fractions at 
the Tevaron and LHC step by step. To further explore the uncertainty from the NRQCD 
factorization scale /xa dependence and different choice of the Xb(3P) feed-down ratios, we 
have performed fits in three schemes by using different NRQCD factorization scale pa and 
Xb{3P) feed-down ratios. 

All the obtained results can explain the transverse momentum distribution of production 
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rate very well just as in our previous work. And for polarization, the results of T(3S') at 
large pt can reproduce the experimental data, which proves our emphasis in our previous 
work that the the Xb(3P) feed-down could be very important for T(3S') polarization. The 
behavior for T(1S', 2S') changed little with or without Xb(3P) feed-down. And the polariza¬ 
tions results can explain the CMS data well, but the distance from CDF data can not be 
smeared. The results for the fraction show different behaviours for each n = 1, 2, 3, 

they are increasing, flat {pt >10 GeV), slightly decreasing respectively as pt increases, and 
all of them £t the experimental data well. We also presented our prediction on the ra¬ 
tio (j[xf, 2 (l-P)]/o'[X 6 i(lP)]; which can reproduce the CMS measurements well, but a little 
underestimated the LHCb data. 

In our study on the uncertainty with three schemes, we hnd that the different choice of 
Xfe(3P) feed-down ratios almost did not modify the hnal results since its effect is almost 
renormalized by the its CO LDME which is the dominant over its color- 

singlet part (the value is hxed). Therefore this study can not distinguish which choice of 
Xfe(3P) feed-down ratios is better, only with further experimental measurement on Xb(3P) 
hadroproduction the feed-down ratio can be hxed. 

Furthermore, we hnd that uncertainty from the NRQCD factorization scale depen¬ 
dence is viewable in our ht and predication. The important fact is that there are sizeable 
diherence for the obtained CO LDMEs in the hts between diherent choice of NRQCD fac¬ 
torization scale although the htted result is almost the same. The diherent CO LDMEs sets 
can bring much diherent prediction at other experiments such as ee, ep colliders. The un¬ 
certainty is due to that the matching between NRQCD and QCD can not be exactly made in 
the hxed-order perturbative calculation, which has also been observed and discussed in pre¬ 


vious work 


u 


35l | . Therefore it must be taken into consideration in the result presentation 


or global ht. 
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